Universal Behavior of Charged Particle Production in Heavy Ion 
Collisions at RHIC Energies 

Peter A. Steinberg 2 * for the PHOBOS collaboration 

B.B.Back 1 , M.D.Baker 2 , D.S.Barton 2 , R.R.Betts 6 , M.Ballintijn 4 , A.A.Bickley 7 , R.Bindel 7 , A.Budzanowski 3 , W.Busza 4 , 

A. Carroll 2 , M.P.Decowski 4 , E.Garcia 6 , N.George 1 ' 2 , K.Gulbrandsen 4 , S.Gushue 2 , C.Halliwell 6 , J.Hamblen 8 , 
G.A.Heintzelman 2 , C.Henderson 4 , D.J.Hofman 6 , R.S.Hollis 6 , R.Holynski 3 , B.Holzman 2 , A.Iordanova 6 , E. Johnson 8 , 
J.L.Kane 4 , J.Katzy 4 - 6 , N.Khan 8 , W.Kucewicz 6 , P.Kulinich 4 , C.M.Kuo 5 , W.T.Lin 5 , S. Manly 8 , D.McLeod 6 , 
J.Michalowski 3 , A.C.Mignerey 7 , R.Nouicer 6 , A.Olszewski 3 , R.Pak 2 , I.C.Park 8 , H.Pernegger 4 , C.Reed 4 , L.P.Remsberg 2 , 
M.Reuter 6 , C.Roland 4 , G.Roland 4 , L.Rosenberg 4 , J.Sagerer 6 , P.Sarin 4 , P.Sawicki 3 , W.Skulski 8 , S.G.Steadman 4 , 
P.Steinberg 2 , G.S.F.Stephans 4 , M.Stodulski 3 , A.Sukhanov 2 , J.-L.Tang 5 , R.Teng 8 , A.Trzupek 3 , C.Vale 4 , 

G.J. van Nieuwenhuizen 4 , R.Verdier 4 , B.Wadsworth 4 , F.L.H. Wolfs 8 , B.Wosiek 3 , K.Wozniak 3 , A.H.Wuosmaa 1 , 

B. Wyslouch 4 

Argonne National Laboratory, Argonnc, IL 60439-4843, USA 2 Brookhaven National Laboratory, Upton, NY 
11973-5000, USA 3 Institute of Nuclear Physics, Krakow, Poland 4 Massachusetts Institute of Technology, 
Cambridge, MA 02139-4307, USA 5 National Central University, Chung-Li, Taiwan 6 University of Illinois at 
Chicago, Chicago, IL 60607-7059, USA 7 University of Maryland, College Park, MD 20742, USA 8 University of 
Rochester, Rochester, NY 14627, USA 



The PHOBOS experiment at RHIC has measured the multiplicity of primary charged particles as a function 
of centrality and pseudorapidity in Au+Au collisions at ^/sjvjv = 19.6, 130 and 200 GeV. Two kinds of universal 
behavior are observed in charged particle production in heavy ion collisions. The first is that forward particle 
production, over a range of energies, follows a universal limiting curve with a non-trivial centrality dependence. 
The second arises from comparisons with pp/pp and e + e~ data. {N c h) / (N pa rt/2) in nuclear collisions at high 
energy scales with y/s in a similar way as N c h in e + e~ collisions and has a very weak centrality dependence. 
This feature may be related to a reduction in the leading particle effect due to the multiple collisions suffered per 
participant in heavy ion collisions. 



1. Introduction 

The PHOBOS experiment has measured 
dN c h/drj and the average multiplicity of charged 
particles (N c h) produced in heavy ion collisions 
for center of mass energies in the nucleon-nucleon 
center of mass system, ^/snn, of 19.6, 130 and 
200 GeV. The data is also binned as a function of 
event centrality (impact parameter) characterized 
by the number of participating nucleons, N part , 
allowing comparisons to elementary systems, like 
pp/pp and e + e~ — * hadrons. 

The PHOBOS multiplicity detector consists of 
several arrays of silicon detectors which cover 
nearly the full solid angle for collision events. The 
event centrality is characterized by the multiplic- 
ity of charged particles measured by two sets of 
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16 paddle counters covering 3 < \r)\ < 4.5. The 
methods used for measuring the multiplicity of 
charged particles as well as for extracting N part 
has been described in more detail in Ref . jj] . 

2. Limiting Behavior in Pseudorapidity 
Distributions 

Fig. | shows dN ch /dr)'/(N part /2) (r)' = r) - 
2/beam) measured at three different RHIC energies 
for peripheral ((N part ) ~ 100) and central events 
((N par t) ~ 355), in the left and right panels, re- 
spectively. These show a clear "limiting behav- 
ior" in the fragmentation region. That is, the 
distributions are independent of beam energy in 
a substantial range in rf . As the beam energy 
increases, dN/drf follows the universal trend un- 
til it reaches 85-90% of it's maximum value at 
midrapidity, at which point it stops following the 
trend. Similar behavior has been observed in el- 
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ementary collisions as well, both in pp collisions 
III and in e + e~ collisions over a large range of 
energies M. 
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Figure 1. dN/drf / (N part /2) for peripheral 
((Npart) ~ 100) and central ((N part ) ~ 355) 
events at three RHIC energies. 

The limiting curve constrains the energy de- 
pendence of the charged particle multiplicity. It 
also varies with centrality in such a way that the 
increases seen at low 7/ (which is midrapidity in rj) 
as Np ar t increases, are accompanied by decreases 
near 7/ ~ (forward rapidities), as seen in Fig. 
[|. It is not clear why this behavior occurs, e.g. 
whether it is from energy conservation or a true 
long-range correlation. 

3. Comparison with Elementary Systems 

The angular distributions of charged parti- 
cles for different strongly-interacting systems are 
shown in Fig. |[ where central Au+Au (divided 
by (Np ar t/2)), pp H and e + e~ (J] data are com- 
pared, all at y/s = 200 GeV. In this figure, the 
Au+Au and pp data are shown as dN c h/dr], while 
the e + e~ data (with cuts applied to reject events 
with large initial-state photon radiation) is shown 
as dN/dyr, the rapidity distribution along the 
event thrust axis, calculated assuming the pion 
mass. The shapes of Au+Au and e + e~ are sim- 
ilar (within 10%) in shape and magnitude, espe- 
cially within |?7| < 4. In the lower panel, we also 



Figure 2. Top: dN^/di]/ (Np ar t/2) for central 
Au+Au collisions at y/SNN = 200 GeV compared 
with pp and e + e~ data. Bottom: Au+Au and pp 
data divided by a fit to the Au+Au data. 

observe that the shapes of the Au+Au and pp dis- 
tributions are also very similar over a large range 
in 77, but the integrals differ by ~ 40%. 

In Fig. |, we compare (N c h)/(N part /2) in 
heavy ion collisions Q to e + e~ and pp/pp data 
over a large range in ^/s |||. It is observed that 
(N c h) I '(Npart /2) lies below pp at low energies, 
passes through the pp data around y/s ~ 10 GeV, 
and then gradually joins with the e + e~ trend 
above CERN SPS energies. These comparisons 
can be seen more clearly by dividing all of the 
data by a fit to the e + e" data [Q. 

The pp/pp data follows the same trend as e + e~ , 
but it can be shown that it matches very well 
if the "effective energy" v^off = is used, 

which accounts for the leading particle effect seen 
in pp collisions Q . Ref. § finds that bulk parti- 
cle production in pp and e + e~ data does not de- 
pend in detail on the collision system but rather 
the energy available for particle production. In 
this scenario, the Au+Au data suggests a sub- 
stantially reduced leading particle effect in central 
collisions of heavy nuclei at high energy. 
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Figure 4. (N ch )/(N part /2) vs. 
v /s^F = 19.6, 130, and 200 GeV. 



Npart for 



Figure 3. Comparison of (N c h) / (N part /2) for 
A+A, pp/pp, and e + e~ data compared with a fit 
to the e + e~ data. 



The alleviation of the leading particle effect 
might not be so surprising in nuclear collisions. 
Each participating nucleoli is typically struck 77 = 



N "" t "i2 > 3 times on average as it passes through 
the oncoming gold nucleus for N part > 65. One 
could speculate that the multiple collisions trans- 
fer much more of the initial longitudinal energy 
into particle production. This naturally leads to 
the scaling of total particle production in heavy 
ion collisions with N part , as seen in Fig. 0, rem- 
iniscent of the "wounded nucleon model" || but 
with the scaling factor determined by e + e~ rather 
than pp. 

In conclusion, PHOBOS has observed two 
kinds of universal behavior. The first is an 
energy-independent, but centrality-dependent, 
universal limiting distribution of charged parti- 
cle production away from midrapidity. The sec- 
ond is that the total charged particle multiplicity 
per participant pair in heavy ion collisions above 
CERN SPS energies scales with sqrt(s) in a sim- 
ilar way as e + e~ collisions. These features may 
offer a new perspective on particle production in 
heavy ion collisions. 
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